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c | ■ Abstract 

P-i! 

Gamma ray burst (GRB) objects are now widely thought to be at 

cosmological distances, and thus represent enormous energy emission. 

^SZ3 ■ Gamma ray spectra extending to GeV energies suggest the possib- 

lity of accompanying neutrino emission, and there are several models 
proposed suggesting the potential detectability of such coincident neu- 
trino bursts. With this in view, we examine possible measurements 
that might be conducted to give experimental data useful for astron- 
omy, for cosmology and also neutrino properties. 

Of interest to astronomy and cosmology, we show how measure- 
ment of neutrino flavor ratios yields information on the nature and rel- 
ative distance of the source. We point out that cosmological time dila- 
tion might be measured for these sources using neutrinos, as has been 
done for photons, and that neutrino oscillation lengths in the range 
of 1 to 10 5 Mpc can be probed with GRB neutrinos. We thus note 
that these sources may make possible the first non-electromagnetic 
measurements of the scale size of the universe. We discuss tests of 
the weak equivalence principle, tests for flavor dependent gravita- 
tional couplings, and tests for long time scale variation of physical 



constants. We also show that a number of new bounds on neutrino 
properties (charge, mass, speed, lifetime) could be facilitated to levels 
well beyond those already inferred from the neutrino observation of 
SN1987A. 

We also examine the implications of these physics opportunities for 
designers of neutrino telescopes. We conclude that detection may be 
possible in planned instruments if the spectra are power law extending 
to the TeV energy region, and if the neutrino fluxes are equal to or 
greater than the gamma ray fluxes. We emphasize the importance 
of low energy detection in future experiments for the tests described 
above. 



1 Introduction 

Recent observations of gamma ray bursts (GRBs) supports the idea that 
these objects may be distributed over cosmological ranges 0. The GRBs may 
be among the most luminous objects in the universe, with peak gamma ray 
luminosities L ~ 10 51 ergs/ sec. There is evidence that the range of intrinsic 
luminosities is very narrow, which might allow distance determinations out 
to red-shifts of 2 or more0. It is plausible that when such large amounts 
of energy are released in objects which are evidently quite compact, pions 
will be produced in hadronic collisions, and that the GRBs may serve as 
standard candles in bursts of neutrinos as well. The (~ 1/E 2 ) power law 
spectrum of gamma rays which is observed to extend at least to several GeV , 
suggests the possibility of particle acceleration and radiation from other than 
electromagnetic interactions of electrons. Paczynski and Xu |§ have recently 
proposed a specific GRB model in which the gamma rays and neutrinos 
arise from decay of pions produced in shock front collisions. In a second 
class of models (e.g. Plaga |1) gamma rays are hypothesized to come from 
superconducting cosmic strings (SCSs); the luminosities are very high and 
one expects neutrino emission as well. 

Models of the first kind[|3| will have the following generic features for the 
neutrinos emitted: the neutrino energies range over MeV to (a few) GeV, or 
perhaps much higher, and since the source is 7r-decay, the flavor content has 
the proportions u e : v^ : v r — 2 : 1 : 0. We also note that due to gamma ray 
absorption within the source, the luminosity in neutrinos can, in principle, 



exceed the observed luminosity in gamma rays, potentially by a large factorQ 
In models of the second kind[|J, the neutrino flavor content is expected to be 
v e : z/ M : u T =1:1:1 and the energies can reach up to 10 TeV . 

In this paper we suppose that the GRBs are indeed cosmological standard 
candles (at least in an ensemble average sense) for neutrinos and gammas, 
in order to maximize the physical inferences to be drawn from coincident 
neutrino and photon detection. We investigate the opportunities for neutrino 
astronomy, for cosmology, and for particle physics. Some of this analysis may 
also be applicable to other astrophysical neutrino sources, such as AGNs 
which are expected to emit neutrinos over a very large energy range; AGNs 
will probably have observational data exceeding what is possible from GRBs, 
but are unlikely to have such short time pulse emission^. We hope that our 
analysis will be useful to designers of neutrino telescopes and inspire them to 
consider GRB neutrino detection as an important goal for future instruments. 

2 Neutrino Mass Mixing and Oscillations 

One of the unique properties of three generations of neutrinos, not shared 
with photons, is that they may carry with them two additional scale lengths 
associated with flavor oscillations. The oscillation lengths are related to the 
neutrino mass differences according to 

; « (km) ^- 5| ^P)i- (1) 

where i and j are flavors and 5rriij is the, currently unknown, mass differ- 
ence between neutrinos of the relevant flavors. The probability of observing 
neutrinos of a given flavor at any distance from their source is a well known 
function of the distance, the oscillation length and mixing parameters, as well 
as the initial conditions. Therefore, if the neutrinos are produced in the same 



1 Recent emphasis has centered on neutron star models H, motivated in part by the 
observation that the gamma-ray energy released by GRBs at cosmological distances is of 
order of a per cent of the binding energy of a neutron star if the emission is isotropic, and 
less if the emission is beamed. In analogy to a supernova one might expect 99 to 99.9% 
of the binding energy is expected to emerge as neutrinos, yielding a neutrino-to-photon 
energy emission ratio of 10 2 to 10 3 . 

2 However, Markarian 421 has been seen by the Whipple Observatory very-high-energy 
7-ray telescope to have a flux increase of a factor of ten on a time scale of one 



flavor ratio in each GRB, then the measured flavor ratios can, in principle, 
determine the relative distances between various sources, at least for the case 
of small mass differences and consequent long oscillation lengths. Absolute 
distances will not be given directly by the measured flavor ratios, because 
the oscillation length itself will be difficult to measure independently. 

However, if GRBs are standard candles whose distances can be calibrated 
(e.g., by the time dilation of their spectra, or by gravity waveforms of inspi- 
raling objects^), then the measured flavor ratios may provide a direct deter- 
mination of very small neutrino mass differences. Competing models place 
the GRBs (i) within the solar neighborhood (D ~ pc), (ii) in the galactic 
halo (I0~ 2 kpc ~ D ~ I0 2 kpc), and (hi) at cosmic distances (D ~ Mpc) as 
we assume herein. Assuming detection of a correlated neutrino with GeV 
energy, the £m 2 's which are probed at each distance scale are ~ 10~ 13 eV 2 , 
~ KT 14 to ~ 10- 18 eV 2 , and ~ 10~ 19 eF 2 , respectively. 

If the neutrino oscillation indications from the atmospheric neutrino ob- 
servations are correct {8m 2 ~ 0.01 eV 2 , and MCC large, ~ 20° — 40° Q), 
then the observed flavor content of the neutrinos from pion-producing GRBs 
should show the same effect, viz. a v^ : v e flavor content of 1.2:1 rather than 
2:1. This applies to models of the Paczynski-Xu kind. Furthermore, because 
the atmospheric solution has a short oscillation length l^ ~ 250p(GeV) km, 
this flavor ratio should be universal for all GRBs independent of their dis- 
tance; individual GRBs might have hot source spots smaller in size than the 
oscillation length I, but individual oscillations will sum to zero in an ensemble 
average. For the models like that of Plaga, based on SCS's, the democratic 
flavor mixture is unchanged by oscillations and remains v^ : v e : v T =1:1: 
1. It may be thus possible to distinguish between these two classes of source 
models. 

If the atmospheric neutrino anomaly is not due to neutrino oscillations, 
then there is a more interesting possibility for GRB neutrinos. According to 
the see-saw mechanism ||, 

m v ~ m 2 D /M 

where mo is the generation-dependent Dirac mass of either the up quarks, 
down quarks, or charged leptons mi. If we use mn ~ mi and M ~ Planck 
Mass, the neutrino masses]] are m Ue ~ 2 ■ lO -17 ^, m v ~ 10 _12 eV r , m Vr ~ 

3 Neutrino mass differences induced by scattering on the inter-galactic medium are 



3 ■ 10- 10 eV, and hence 5m 2 efl ~ 10~ 24 eV 2 and 5m\ r ~ 10- 19 eV 2 . 

In this see-saw case, the oscillation lengths are in the cosmic range 1 
to 10 5 Mpc for neutrino energies in the GeV rangeQ. There is a maximum 
distance over which a mixture of mass eigenstates can be expected to remain 



coherent [|TT|1 . However, for the tiny mass differences here, the coherence 
length appears to exceed 10 5 Mpc. 

The longer oscillation length possibilities exceed the size of the visible 
universe ~ Hq = 3000/i _1 Mpc, where Hq = lQOh km/sec/Mpc is the 
present value of the Hubble parameter. If / exceeds Hq 1 , then the neutrinos 
do not oscillate, and the flavor ratios observed at earth are just those estab- 
lished at emission. On the other hand, for I less than ~ 10 4 Mpc, the flavor 
content of the neutrinos from pion-producing GRBs will vary strongly with 
distance; and again, for Plaga-like models the flux will remain universal with 
no dependence on distance. 

Alternatively, if mass differences are large (of order 0.01 eV 2 or larger), 
then the separation of the mass eigenstates in time offers another handle on 
5m 2 . At a fixed energy E, the arrival times of the z/'s are separated by 

,(L/100^)(faVlO-'eV) 

(E/100 MeV) 2 

assuming much smaller time difference at the source. For low energies and 
large 5m 2 , 5t can range from 10~ 3 sec to several seconds. With a spectrum 
of energies, the neutrino burst would be spread more than the accompanying 
gamma ray burst. 



3 Other Neutrino Properties 

As with the observations of Supernova 1987A, the distance, short emission 
time and trajectory through varying gravitational fields, leads to the poten- 
tial for some fundamental tests of neutrino properties, which are not possible 
in terrestrial laboratories |12[. In the following we enumerate some of the 
possibilities. 



expected to be smaller than these numbers) 10 . 

4 For most purposes in this paper, it is sufficient to equate time and distance as if the 
universe were static. The generalization to an expanding cosmology is reserved for §4. 



3.1 Neutrino Lifetime 

From the detection of neutrinos arriving form the 50kpc distant SN1987A it 
was possible to place a limit on the (laboratory frame) lifetime of the i> e of 
r(i7 e ) > 5 • 10 12 sec for E v of the order of 10 MeV. In the same way, ob- 
servation of z/'s from GRBs would place bounds on t{v) > 10 17 (D / Gpc) sec, 
which for a 10 Mpc distant source is some 200 times stronger than the bounds 
on z7 e from SN1987A. Of course this is really a bound on the lifetime of the 
dominant mass eigenstate. 

3.2 Neutrino Electric Charge 

If neutrinos carry even a tiny electric charge (as favored in some theoretical 
models JT3J), then the passage through magnetic field regions enroute to the 
earth from a distant source creates an additional dispersion in arrival time, 
St. From the observed upper bound on any additional dispersion beyond the 
gamma pulse time distribution one can then derive a limit on the neutrino 
charge from the Barbillieni-Cocconi formula [|14||: 




-^ (5E/E)-^ 
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where St is the burst duration, D is the flight distance, B is the root- 
mean-squared average magnetic field experienced by the neutrinos enroute, 
and SE/E is the relative spread in energies. For St ~ 10~ 3 sec, D ~ 10 Mpc, 
< E >~ 1 GeV, B ~ 10~ 12 Tesla and SE/E ~ 1, one has Q u /\e\ < 10" 27 . 
Hence it may be possible to improve considerably on the SN1987A limit of 
Q Ue < 10~ 14 |e|, and the somewhat better laboratory limit of Q Ve < 10~ 19 |e|. 
Moreover, one could place the first limit on v^, and possibly on v T . 

3.3 Neutrino Speed 

To the extent that the 7-ray and neutrino pulses coincide, new limits on 
neutrino speed relative to the speed of light may be placed. If St = t u — t 1 , 
then 

1 - v„/c < St/D. 



For St as large as 1 second this can place upper bounds of the order 
of 10~ 15 (lOMpc/ D) on 1 — (3 V , to be compared with 10~ 9 (for u e ) from 

SN1987AI 



3.4 Equivalence principle 

The neutrinos from SN1987A were used to establish limits on parameters 
in various theories of gravitation and to test whether the Weak Equivalence 
Principle (WEP) is symmetric with respect to bosons and fermions, and 
with respect to matter and anti-matter. Specifically, the Shapiro delays for 
gammas, neutrinos, and anti-neutrinos passing the Galactic nucleus were 
compared and found to be the same, within errors |TB|, IT7J. Also, SN1987A 



provided tests for the presence of proposed new forces in nature and tests 
of the equality of parameters when applied to photons, neutrinos, or anti- 
neutrinos [pj], |I8|| . The same methods would apply to neutrinos from GRBs. 
However, under our assumptions, the distances would be greater and the 
impact parameter with the Galactic nucleus would vary from event to event, 
offering much improved sensitivity and a new distance scale, as well as a 
large increase in statistics. 

To test the possibility that the neutrino and antineutrino gravity cou- 
plings differ requires tagging of neutrino and antineutrino events separately. 
Distinguishing between v e and v e interactions at low energy may be possi- 
ble because of the charged current capture by protons. At higher energies 
distinguishing between v^ and v^ interations is in principle possible if one 
employed a magnetic field, or if one could adequately detect the muon cap- 
ture by nuclei; but in practice it cannot be done in instruments proposed at 
this time. A v — v separation would also be useful for studying possible CP 
violation in neutrino oscillations |BJ. 

A flavor violating gravitational coupling has been proposed as a pos- 
sible mechanism for accounting for atmospheric neutrino as well as solar 
anomalies [pOJI - It is remarkable that one choice of parameters (sin 2 26G ~ 1, 
5f ~ O(10~ 15 )) can account for both problems. The transition probability 
for v^ <^> v e is sin 2 (29G)sin 2 (^LE(j)(L)Sf), where is the gravitational po- 
tential and Sf is a measure of the degree of violation of the WEP. Without a 
knowledge of <j) along the path it is not possible to calculate the net effect on 
the mixing, but it is very likely similar to the expectation in the oscillation 
case (i.e., v e : v^ becomes 1.2 : 1 from 2:1). 



4 Cosmology 



Norris, et al. [pl| , have analyzed the experimental data assuming that the 
GRBs are standard candles in gamma rays and have found evidence for a 
cosmological time dilation. If the GRBs are neutrino sources, then the same 
analysis can provide, for the first time, non-electromagnetic evidence for the 
expansion of the universe^]. 

The prevailing view, that the expansion of the universe is of a universal 
nature, leads to the expectation that the gamma ray and neutrino time dila- 
tions will be identical. Unfortunately, a nearly identical dilation may occur 
if the cause is evolutionary effects. This is because the charged and neutral 
pions of the pion-pro duct ion model, and the emitted photons and neutrinos 
of the SCS model arise from a common mechanism. Still, it is possible that 
the neutrino and photon opacities evolve differently, in which case studies of 
the differing dilations may yield information on cosmic evolution. 

An important probe of the universe at an early time is the oscillation 
phase itself, 0y. In Minkowski space, the quantum mechanical phase is just 
Et. However, it is a bit more complicated in the expanding universe. In the 
adiabatic approximation, the phase is 

r T (It' 

*=*°Lwj E - (2) 

We have assumed a time-dependent Robert son- Walker metric, with scale 
factor .R(t); r is the lookback time to the source emission, and R(0) = Rq. 
The red-shift factor R /R(t) accounts for the time dilation of our clocks com- 
pared to early universe clocks. Taking the difference of two mass-eigenstate 
phases, expanding Ei ~ p + mf/2p, and red-shifting the momentum to its 
present-time observed value via p(r) = (Ro/R(r))p , one obtains from Eq. 
(PD the very simple result for the oscillation phase: 
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The blue-shift of the inverse momentum exactly cancels the red-shift 
from time dilation. Thus, the correct generalization of the oscillation phase 



5 We note that such an analysis would exclude any cosmological model which attempts 
to explain the observed red-shifts by "tired photons" ; earlier arguments against the tired 
photon hypothesis were given by Zeldovich, et a/.p2|. 
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from Minkowski space to an expanding cosmology is obtained by replacing 
laboratory time with cosmic lookback time, or equivalently, replacing Uj in 
Eq. ([]]) with CTij. 

This deceptively simple result is rich in cosmological information. In 
particular, if 5m 2 were a priori known, then a measurement of the flavor 
ratio and E v would directly yield r. This is analogous to obtaining the 
cosmic red-shift z directly from a measurement of a photon's energy, when 
the unshifted spectral line of the photon source is a priori known. Further- 
more, r contains as much information as is conveyed by z. In fact, r, z, 
and the distance D are linearly related for small z, r, D, and nonlinearly re- 
lated for large z,t,D. For example, to first nonlinear order, some Taylor 
series expansions relating these three variables are[E3 
y)(H q t) 2 + . . ., and its inverse relation r(z) = Hq 1 

and Dl = Hq 1 z + |(1 — qo)z 2 + . . . . Here, D L is the "luminosity" distance 
defined as D\ = LjAmT, where £ is the absolute luminosity at the source 
(energy /time), and T is the fluence measured at earth (energy /time/area). 
Results for other distance definitions are similar; e. g. the "proper" distance 
(RoX comoving cordinate distance), is given as Dp = Hq 1 z — \{1 + qo)z 2 + 
In all of these relations, go is the present value of the deceleration parameter. 

The value of go is unknown; since the GRBs seem to exist at cosmic 
distances, it is possible that oscillation measurements with GRB neutrinos 
may shed some "neutrino light" on this important parameter. An idependent 
measurement of even two of the three variables z, r, and D would potentially 
determine go an d test cosmological models, because of the nonlinear relations. 
The series expansions make it clear that the linear Hubble relation fails by 
a fractional amount z in red-shift, Hqt in lookback time, and D/Hq 1 in 
distance. Interpreting the burst dilation of fainter GRBs as due to time 
dilation leads to the estimate z ~ 1 to 2 for these GRBs, so there is indeed 
hope that oscillation measurements may yield a value for g . We have seen 
that for small neutrino masses, it may be possible to measure the phase of 
oscillations from pion-producing GRBs. Of course, one would require a much 
better understanding of the structure and mechanisms of GRBs than we have 
today, in order to draw useful conclusions. 

Let us assume for the moment that such measurements can be made with 
enough precision to yield values for rSm 2 -. If an independent measurement 
of z or D is available, then a single GRB would suffice to fix 5m 2 , and a 



second GRB measurement would yield go- in this idealized situation, we 
may also inquire about the nature of higher order terms in the nonlinear 
expansions relating z,t, and D. In fact, given a cosmological model, the 
nonlinear relations among z, r, and D are exactly calculable. For example, 
r = f*° T) dR/R = fKfodlnR/H will yield r(l + z = Rq/R(t)) once R(t) 

or equivalently, H = R/R are determined as a function of R (or z) from 
the Friedmann equation. Ignoring the radiation energy density of the re- 
cent universe compared to the matter density, the Friedmann equation reads 
H (z) 2 = Hi [(1 + z) 2 (l + zQ ) — z(2 + z)Qa\, where Q is the present mat- 
ter density compared to the critical value p c = 3Hl/8nG, and fl A = A/3Hq; 
A is the comsological constant. (The Friedmann universe is open or closed 
according to whether Q + Q A is less than, or greater than, unity.) The 
Friedmann equation may be manipulated to yield 
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(4) 



where the curvature term Q k = ^o + ^a — 1 is a constrained by the 
Friedmann equation itself. A thorough set of numerical solutions to this 



equation may be found in ref.[2^]. With £7a omitted, the integral is easily 
solved analytically. The form of the solution for t(z) depends on whether 
the universe is closed (f2 > 1), critical (Q = 1) or open (Q < 1). For the 
critical case, motivated by inflation, the result is 
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(5) 



For small z, the integral is just z, and the linear Hubble relation H t ~ z 
results. A related calculation yields the simple and useful relation between 
the red-shift and the luminosity distance [p5l 
(2Ho 1 /n 2 ) \zQ + (fio - 2)(v^no" 



-hi-: 

2H, 1 



valid for ^a <t. ^o : Dl = 
Then for Q = 1 as suggested 



l-Vz+1 



Given a cosmo- 



in inflationary cosmologies, Dl 

logical model, measurements of a second GRB would potentially validate or 

invalidate the model, by fitting or not fitting the nonlinear r(z) or t(D) 

relation. 

In static Minkowski space, a neutrino source is most useful for oscilla- 
tion studies if its distance is comparable to the oscillation length; a shorter 
distance does not provide sufficient path length for oscillations to develop, 



10 



and in a longer distance the information is effectively averaged over many 
oscillations. However, in an expanding metric, distance must be carefully 
defined, and the situation is different. We may calculate the lookback times 
and luminosity distances to typical GRBs using the t(z) or D L (z) formulae 
just given. The results are t(z = 1) = OASHq 1 and r{z = 2) = 0.54i7j" 1 , 
and D L (z = 1) = 1.18i7 _1 and D L (z = 2) = 2.5AH \ Note that for this 
matter-dominated, critical-density example, D L (z) is unbounded as z in- 
creases, whereas the lookback time r(z) has approaced \Hq for large z. 
The lookback time is the relevant variable for neutrino oscillation, since it is 
directly proportional to the oscillation phase, according to Eq. (|3|). The fact 
that t has become asymptotic for the red-shift values typical of GRBs has 
an important implication: the value of r in the oscillation phase is nearly 
t(z = oo), which in any cosmological model is a known fraction of Hq 1 . This 
fact means that the uncertainty in r is dominated by the uncertainty in H , 
which is only a factor of two. Thus, a single measurement of the oscillation 
phase and neutrino energy may permit a determination of 5m 2 to a factor of 
two! 

Comparing an exact cosmological solution for z, r, or D to the appro- 
priate Taylor series approximation relates each term in the series to more 
fundamental quantities. The simplest relation is between the parameter q 
multiplying the quadratic term, and the parameters Qo and ^a in the Fried- 



mann equation. The relation is|24| go — ^o/2 — ^a- Thus, a neutrino oscil- 
lation determination of go would establish a fundamental constraint between 
the two parameters of the Friedmann universe. 

Studying cosmology by measuring neutrino flavor ratios has one tremen- 
dous advantage over other methods, namely that flavor ratios should be 
independent of any evolutionary effects in the ensemble of sources. Abso- 
lute neutrino luminosities may evolve, but the initial neutrino flavor ratios 
are fixed by microphysics; it is hard to imagine that these ratios will change 
with cosmic history. This is in sharp contrast to the use of distant "can- 
dle" luminosities to infer deviations from the linear Hubble Law, which may 
be due to the deceleration of the universe or to evolutionary effects in the 
candles. 



ff 



5 Physical Constants - Time Dependence 

There is a long tradition in physics of asking whether various constants are, 
in fact, constant over cosmological time. The best known of these questions 



was raised by Dirac. In 1981, Barrow [26] reviewed these ideas. Among his 



conclusions is the idea that only dimensionless constants can have a meaning- 
ful time dependence. In fact there are some new strong limits upon the time 
variation of the fine structure constant and the electron to proton mass ratio, 
and upon posible variations across causally disconnected regions of space [f2"7j]. 

However, any possible time variation of the dimensionless parameters 
upon which neutrino oscillations depend (mass ratios and mixing angles, or 
equivalently, the Yukawa couplings at the origin of fermion mass generation) 
are without constraint. If the distances to the GRBs can be measured by 
independent means, and the neutrino mass differences turn out to be very 
small, then a time dependence of the mixing angles over cosmological times 
may be detectable as deviations from the expected flavor ratios. 

It has been speculated that dimensionful comological parameters may 
have a time-dependence. After all, the Hubble parameter itself has a com- 
plicated dependence on time, varying inversely during power law expansion, 
and exponentially during inflationary expansion. The cosmological constant 
A(t) (or equivalently, the energy density of the vacuum), has been hypothe- 



sized to relax to zero asymptotically with time |pq| . And the Hubble param- 
eter and Newton's gravitational constant have been hypothesized to consist 
of two terms each, one standard and one oscillating periodically in time [^] . 
Motivation for the relaxation of A is that the present value of A/8irG is 
known to be less than 10~ 47 Gel/ 4 [^3|, whereas there is no good theoretical 



understanding of why this should be so small. Motivation for the periodic 
oscillation in H and G is that periodic modulation offers an explanation 
for the controversial observation pU[ of a 128/h Mpc quasi-period in deep 
"pencil beam" surveys of galaxy positions. If A(r),G(r), and/or H(t) are 
time-dependent, then the Friedmann equation is modified. As discussed in 
the previous section, the oscillation phases of neutrinos emitted at large look- 
back time r are sensitive to the parameters in the (now modified) Friedmann 
equation. Thus, measured flavor ratios can offer information on the time 
dependences of these cosmic parameters. 
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6 Implications for Neutrino Telescopes 

We discuss some implications of our analysis for designers of neutrino tele- 
scopes. The spectrum of gammas from GRBs has been seen out to energies of 
a few GeV, so lacking a detailed GRB model one would do well to focus upon 
neutrino detection in a similar energy range. There are already weak limits 
on the neutrino flux associated with GRBs from the 1MB experiment [BT]] and 
others. The largest deep mine experiment yet planned, the SuperKamiokande 
detector with a 50 kiloton sensitive volume (scheduled for operation in 1996) 
will have about ten times the sensitivity for events with energies between 
5 MeV and a few GeV as had previous instruments. 

Much further progress for sensitivity down to the few MeV region is not 
presently on the horizon. One possibility would be parasitic use of a megaton 
size detector contructed for observation of supernova neutrinos out to a few 
Mpc. The capability to sense neutrino versus antineutrino interactions via 
interaction characteristics, muon absorption, or magnetic fields to distinguish 
charge, would provide a powerful tool for many of the tests described above. 

However, if the GRB neutrino spectrum extends to energies of many 
GeV or even TeV, then there is more hope for the near future. This 
is because the detectability of signals rises strongly with energy. For ex- 
ample, the AMANDA, Baikal, DUMAND and NESTOR ice/water instru- 
ments now under construction have effective volumes for 100 GeV neutri- 
nos of order 10 6 tonnes, up by two orders of magnitude from underground 



instruments [p^l . Observation of the ratio of muon charged current events to 
muonless events would potentially allow for discimination between the puta- 
tive neutrino-flavor democratic source and the pion source with ir — > fj, — > e 
decay and v^/v e = 2. 

We can make a rough estimate of counting rates for neutrinos by assuming 
a spectral shape which we take to be l/E 2 , and a gamma flux which we 
take to be 1 7 /cm 2 /burst with energy greater than 1 MeV. The ratio of 
neutrinos to gamma rays, labeled as n, could be zero in the case of a purely 
electromagnetic origin of the gammas, in which case most of the foregoing is 
irrelevant except for the important constraint upon the source model. For 
the situation of particle acceleration with power law spectra, as discussed by 
Paczynski, n ~ 1. n ~ 1 holds if the source is not heavily shielded. Even 



better for our purposes, in situations similar to that expected near AGNpB" 



the attenuation of the gamma rays can be severe, while the neutrinos flow 
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freely and so 77 >> 1, possibly even as large as the 10 3 energy emission ratio 
expected in supernovae. 

We take a simplistic model of a standard 10, 000 m 2 effective area (for 
muon detection) instrument (e.g., AMANDA, Baikal, DUMAND or NESTOR). 
In underwater (or ice) experiments the area grows with energy somewhat, 
but for simplicity we make the conservative assumption that this area in- 
dependent of energy (as in mine detectors). We take the effective detector 
area for neutrinos then as the the muon range, times the effective area for 
muons, times the density of the medium, times Avogadro's Number, times 
the neutrino-nucleon cross section. This turns out to be 90 cm 2 for the nom- 
inal detector at 1 TeV, and it scales roughly as E 2 from cross section and 
range, for energies from 1 GeV to 10 TeV. 

A given detector will have some threshold detection energy, which in 
practice is not a step function, though for simplicity we take it to be so, at 
20 GeV. We also take an arbitrary maximum neutrino energy of 1 TeV. 
With these assumptions we find for the expected number of neutrino inter- 
actions per burst which are below the neutrino detector horizon, the value 
9 ■ 10~ 5 7]/burst in muon neutrinos. If we take the rate of GRB (as now de- 
tected) as about 1/day, then the total expected number of correlated muons 
for this standard neutrino detector is 1.5 • 10~ 2 r]/year. Since the input as- 
sumptions are certainly imprecise to a factor of ten, this could easily be well 
detectable or beyond experimental reach. 

If rj = 1, it is easy to see why existing underground instruments have not 
seen such correlations as yet, despite lower thresholds. The 1MB detector had 
400 m 2 area, 25 times less than our assumed instrument. In fact, one can 
turn this around and ask what limit on 7] is implied by the non-detections 
in present underground instruments, under our assumptions of spectrum. In 
Figure [I], we show the combinations of maximum GRB neutrino energy and 
tj which would lead to one event detected per year in 1MB (400 m 2 ), a next 
generation instrument (10 4 m 2 ), and a hypothetical km 3 detector (KM3). 
One sees that for a TeV GRB neutrino cutoff energy, the 1MB limit on the 
neutrino to gamma ratio is of the order of a few thousand. 

Another approach is to consider the brightest GRBs, and ask whether 
or not multiple events from a single source are possible. If we assume that 
the distribution of GRB gamma fluxes are dominated by spatial distribution, 
then 0.1% of the GRBs will be at 0.1 of the distance of the typical burst, and 
offer 100 times the neutrino flux at earth; this is the famous "3/2 law" . Thus 
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in one year of operation during which there will be about 1000 GRBs, one 
might find a GRB with 9 • 10 _3 ?7 muons. One can conclude that detection of 
multiple muons per GRB is likely in the standard next generation instrument 
only if i] 3> 1 . This is illustrated in Figure |2|, where the diagonal lines indicate 
the T) and maximum neutrino energy combinations needed to see an event 
with ten muons once per two years in the various classes of detectors (and 
roughly one in three such GRBs would have a coincident GRO observation). 

Better possibilities exist with a km 3 scale instrument. Even though cur- 
rent design discussions center around optimization for detection of TeV to 
PeV neutrinos from AGN's, neutrinos from GRB's (as well as other oppor- 
tunities such as dark matter searching and atmospheric neutrino oscillation 
measurement) argue in favor of reducing the threshold to as low a value as 
feasible, say of the order of 10 GeV. Given the lack of understanding of 
the process which generates the GRB's we cannot do much in terms of a 
priori optimization of detector design. Good timing is obviously important, 
but high angular resolution is not as demanding as for neutrino point source 
searches. 

If we take the energy threshold to again be 20 GeV for the (km 3 ) neutrino 
detector with 10 6 m 2 effective muon area, then under the same assumptions 
we employed for the standard next generation instrument we get an expected 
event number of 9 • 10~ 3 r]/burst, or about l.brj/year GRO coincident detec- 
tions. More encouraging yet is the rate of multiples: the brightest 0.1% of 
GRBs might produce the spectacular signature of 90/x's once in two years if 
tj should be 100 and the maximum energy 1 TeV, in which case one could 
begin the studies outlined above. 

7 Conclusion 

Any detection of GRBs in neutrinos would have great significance in under- 
standing these enigmatic objects. What we advertise herein is that moreover 
such a detection can lead to fundamental exploration of neutrino physics, 
astrophysics, and maybe even cosmology. This exploration is not possible by 
any other means we know. 

The implications for the telescope designers then are fairly obvious: make 
instruments with as low an energy threshold as practical, and allow for up- 
grades (in terms of energy sensitivity and capability to resolve neutrino fla- 
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vors) to follow the path of discovery. 
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Figure 1: One event per year in coincidence with GRO for detectors of 400 m 2 
(1MB), 10,000 m 2 (BAND), and 10 6 m 2 (KM3), for various hypothetical 
neutrino to gamma ratios (if) and GRB maximum neutrino energy. The 
upper right hand corner region is already ruled out by 1MB data. 
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Figure 2: 10 muons in one burst for brightest GRBs, seen in detectors of 
400 m 2 (1MB), 10,000 m 2 (BAND), and 10 6 m 2 (KM3), for various hypo- 
thetical neutrino to gamma ratios (77) and GRB maximum neutrino energy. 
The upper right hand corner region is already ruled out by 1MB data. The 
rate may be 0.5/year. 
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